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Abstract 

Recently, we developed a method for calculating the lifetime of the 
particle in the special situation where there is no potential barrier, 
as a first step in our efforts to understand the quantum-mechanics 
of magnetic traps. The toy model that was used in this study was 
physically unrealistic because the magnetic field did not obey Laplace's 
equation. Here, we study, both classically and quantum-mechanically, 
the problem of a neutral particle with spin S, mass m and magnetic 
moment fi, moving in two-dimensions in an inhomogeneous physically 
realistic magnetic field given by 

B = B' ± (xx - yy) + B z. 
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We identify 




which is the ratio between the precessional frequency of the particle 
and its vibration frequency, as the relevant parameter of the problem. 

Classically, we find that when \i is antiparallel to B, the particle is 
trapped provided that K < y^4/27. We also find that viscous friction, 
be it translational or precessional, destabilizes the system. 

Quantum- mechanically, we study the problem of spin S = h/2 
particle in the same field. Treating K as a small parameter for the 
perturbation from the adiabatic Hamiltonian, we find that the lifetime 
T esc of the particle in its trapped ground-state is 



esc ~ 128^2 exp 



2 

K 



where T vi ^ = 27ry tuBq/ ' \i (B' ± ) 2 is the classical period of the particle 
when placed in the adiabatic potential V = /J, |B | . 



1 Introduction. 

1.1 Magnetic traps for neutral particles. 

Recently there has been rapid progress in techniques for trapping samples 
of neutral atoms at elevated densities and extremely low temperatures. The 
development of magnetic and optical traps for atoms has proceeded in parallel 
in recent years, in order to attain higher densities and lower temperatures 
|]> @> 5 IE HI- We should note here that neutral traps have been around 
much longer than their realizations for neutral atoms might suggest, and the 
seminal papers for neutral trapping as applied to neutrons and plasmas date 
from the sixties and seventies. Many of these papers are referenced by the 
authors of Refs. |T|, @, [| • In this paper we concentrate on the study of magnetic 
traps. Such traps exploit the interaction of the magnetic moment of the atom 
with the inhomogeneous magnetic field to provide spatial confinement. 

Microscopic particles are not the only candidates for magnetic traps. In 
fact, a vivid demonstration of trapping large scale objects is the hovering 
magnetic top|], |7], ||. This ingenues magnetic device, which hovers in 
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mid-air for about 2 minutes, has been studied recently by several authors 



10, 11, 12, 13, 14 



1.2 Qualitative description. 

The physical mechanism underlying the operation of magnetic traps is the 
adiabatic principle. The common way to describe their operation is in terms 
of classical mechanics: As the particle is released into the trap, its magnetic 
moment points antiparallel to the direction of the magnetic field. While 
inside the trap, the particle experiences lateral oscillations uo,^ which are 
slow compared to its precession u prec . Under this condition the spin of the 
particle may be considered as experiencing a slowly rotating magnetic field. 
Thus, the spin precesses around the local direction of the magnetic field 
B (adiabatic approximation) and, on the average, its magnetic moment fi 
points antiparallel to the local magnetic field lines. Hence, the magnetic 
energy, which is normally given by — fi • B, is now given (for small precession 
angle) by \i |B|. Thus, the overall effective potential seen by the particle is 

K//^|B|. (1) 

In the adiabatic approximation, the spin degree of freedom is rigidly coupled 
to the translational degrees of freedom, and is already incorporated in Eq.(|l]). 
Thus, under the adiabatic approximation, the particle may be considered as 
having only translational degrees of freedom. When the strength of the mag- 
netic field possesses a minimum, the effective potential becomes attractive 
near that minimum and the whole apparatus acts as a trap. 

As mentioned above, the adiabatic approximation holds whenever uj prec ^> 
u V ib- As Up rec is inversely proportional to the spin, this inequality can be 
satisfied provided that the spin of the particle is small enough. If, on the 
other hand, the spin of the particle is too large, it cannot respond fast enough 
to the changes of the direction of the magnetic field. In this limit u prec -C oo^, 
the spin has to be considered as fixed in space and, according to Earnshaw's 
theorem[n|, becomes unstable against translations. Note also that uj prec is 



proportional to the field |B|. To prevent uj prec of becoming too small, resulting 
in spin-flips (Majorana transitions), most magneto-static traps include a bias 
field, so that the effective potential V e ff possesses a nonvanishing minimum. 
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1.3 The purpose and structure of this paper. 

The discussion of magnetic traps in the literature is, almost entirely, done 
in terms of classical mechanics. In microscopic systems, however, quantum 
effects become dominant, and in these cases quantum mechanics is suited for 



the description of the trap (T^|. An even more interesting issue is the un- 
derstanding of how the classical and quantum descriptions of a given system 
are related. 

As a first step in our efforts to understand the quantum-mechanics of 
magnetic traps, we recently developed a method for calculating the lifetime 
of the particle in the special situation where there is no potential barrier |T7| . 
The toy model that was used in this study consisted of a particle with spin, 
having only a single translational degree of freedom, in the presence of a 
ID inhomogeneous magnetic field. We found that the trapped state of the 
particle decays with a lifetime given by ~ 1/ (y/Kuvibj exp (2/K) where 

K = u V ib / u prec . The field that was used in this model was not divergenceless, 
and in this sense, the model is unrealistic. The next step, presented in this 
paper, is to study, both classically and quantum-mechanically, the case of a 
particle with spin, having two translational degree of freedom, in the presence 
of a physically realistic (i.e. divergenceless) inhomogeneous magnetic field. 



This model is reminiscent of a Ioffe-Pritchard trapg, |18| , but without the 
axial translational degree of freedom. We neglect the effect of interactions 
between the particles in the trap and so we analyze the dynamics of a single 
particle inside the trap. 

The structure of this paper is as follows: In Sec.(0) we start by defin- 
ing the system we study, together with useful parameters that will be used 
throughout this paper. Next, we carry out a classical analysis of the problem 
in Sec. (|3|). Here, we find two stationary solutions for the particle inside the 
trap. One of them corresponds to a state whose spin is parallel to the direc- 
tion of the magnetic field whereas the other one corresponds to a state whose 
spin is antiparallel to that direction. When considering the dynamical sta- 
bility of these solutions, we find that only the antiparallel stationary solution 
is stable. We also study the same problem but with viscous friction added, 
and arrive at the result that friction destabilizes the system. In Sec.(f|) we 
reconsider the problem, from a quantum-mechanical point of view. Here, we 
also find states that refer to parallel and antiparallel orientations of the spin, 
one of them being bound while the other one unbounded. In this case, how- 
ever, these two states are coupled due to the inhomogeneity of the field and 
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we move on to calculate the lifetime of the bound state. Finally, in Sec.(^) 
we compare the results of the classical analysis with these of the quantum 
analysis and comment on their implications to practical magnetic traps. 



2 Description of the problem. 

We consider a particle of mass m, magnetic moment \i and intrinsic spin S 
(aligned with //) moving in an inhomogeneous magnetic field B given by 

B =B z+B' ± (xZ-yy) . (2) 

This field possesses a nonzero minimum of amplitude at the origin, which is 
the essential part of the trap. The Hamiltonian for this system is 

2 

H=2--n-B (3) 

where p is the momentum of the particle. 

The Hamiltonian is invariant under a group of operations consisting of a 
rotation of position space about the z-axis by an arbitrary angle 7 combined 
with a rotation of spin space about the S^-axis by the opposite angle —7. 
Since the generators of these two rotations are the ^-components of orbital 
angular momentum L z = xp y — yp x and of spin angular momentum S z , 
respectively, this symmetry gives rise to a constant of motion, 

A = L z — S z = const. (4) 

Since the magnetic field B does not depend on z, the motion along the in- 
direction is trivial. Therefore, we restrict ourselves to studying the motion 
in the (x, y)-plane. 

We define u prec as the precessional frequency of the particle when it is at 
the origin (x = 0, y = 0). Since at that point the magnetic field is B =B z 
we find that 

^prec = g ■ {£>) 

Next, we define u v n, as the small-amplitude vibrational frequency of the par- 
ticle when it is placed in the adiabatic potential field given by 

V(x) = ^ \B(x) I = ^B h + \ (^j 2 (x 2 + y 2 ))+0 (x\ x 2 y 2 , 2/ 4 ) ■ 
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For this potential we have 



d 2 V {B' L f 

@x=0 -°0 



dx 2 



and therefore 




^ E V- = ^- ( 6 ) 



We also define the ratio between and ca 



preci 



K _ u vib _ / S 2 (B' ± ) 2 

Uprec V /imBo 

This will be our 'measure of adiabaticity'. It is clear that as K becomes 
smaller and smaller, the adiabatic approximation becomes more and more 
accurate. Note that when the bias field B vanishes, K becomes infinite, 
and the adiabatic approximation fails. We will later show that, under this 
condition, the system become unstable against spin flips, which is in agree- 
ment with our discussion at the beginning. This shows that the introduction 
of the bias field B , is essential to the operation of the trap with regard to 
spin-flips. Note also that K is the only possibility to form a non-dimensional 
quantity (up to an arbitrary power) out of the parameters of the system. 
The value of K therefore, completely determines the behavior of the system. 



3 Classical analysis. 

3.1 The stationary solutions. 

We denote by n a unit vector in the direction of the spin (and the magnetic 
moment). Thus, the equations of motion for the center of mass of the particle 

are 

d 2 y d 
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and the evolution of its spin is determined by 



S^ = fihx B. (9) 

It is straightforward to check that the quantity A = L z — S z is indeed con- 
served. 

The two equilibria solutions to Eqs.(H) and are 

n(t) = Ti (10) 

with 

x{t) = 
y(t) = 

representing a motionless particle at the origin with its magnetic moment 
(and spin) pointing antiparallel (n(t) = — z) to the direction of the field at 
that point and a similar solution but with the magnetic moment pointing 
parallel to the direction of the field (n(t) = +z). 



3.2 Stability of the solutions. 

To check the stability of these solutions we now add first-order perturbations. 
We set 

n(t) = Tz+e x (t)±+e y (t)y (11) 
x{t) = + 5x(t) 
y(t) = + 6y(t), 

(note that, to first order, the perturbation 5h = e a; (t)x+ej / (t)y is taken to be 
orthogonal to the value of n for the stationary solution n = =Fz, since n is 
a unit vector) substitute these in Eqs.(||) and (||), and retain only first-order 
terms. We find that the resulting equations for 5x(t), Sy(t), e x {t) and e y {i) 
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are 



d 2 5x _ pB' ± 

at m 
d 2 5y pB' A 



dt 2 ™ -v 



rn 



The normal modes of the system transform as the irreducible representa- 
tions of the symmetry group. The 4-dimensional linear space spanned by the 
deviations (Sx, 6y, e x , e y ) from the stationary state carries the irreducible rep- 
resentations T + with characters e - * 7 and T_ with characters e + * 7 , and may 
thus be decomposed into the two 2-dimensional invariant subspaces trans- 
forming as T + and T_, respectively. These subspaces are spanned by the 
circular position coordinates and precessional spin coordinates 

r+ : (p + = 5x + i5y, e_ = e x - ie y ); (13) 
r_ : (p_ = 8x - iSy,e + = e x + ie y ). (14) 

Thus, the normal modes consist of a circular motion in the (x, y)-plane cou- 
pled to a precession of the spin vector in the opposite sense. 

Indeed, after introducing the (p±, e T )-coordinates into Eqs.(|i~2D, this set 
of four equations decomposes into one pair of equations for (p + ,e_) and 
another pair for (p_,e+). We now look for oscillatory (stable) solutions of 
these equations and set 

P± = Pi.oe-*"*, e± = e ±j0 e-^. (15) 
This yields the algebraic equations 

u 2 mb B /B' ± \ ( p +fi \ _ / 



1 1 a\ 

±i<jj prec B' ± /B i(uj + uj prec ) J V e_ J V 1 ' ' 



Tiu prec B' ± / B 



l[UJ 



prec ) 



P-,0 
e +,0 








(17) 
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These equations have non-trivial solutions whenever the determinant of either 
of the two matrices vanishes. This yields the secular equations 



T + : k(^-Y+(^-)\i = 0, (18) 
r_: K (^)'-(^)'±l=0, (19) 

\^vib/ \^vib/ 

which determine the eigenfrequencies u of the various modes. Since the sys- 
tem has three degrees of freedom, we expect to have three normal modes. 
Indeed, when uj is a solution of the first equation, then — uj is a solution 
of the second equation. We define the mode frequencies in Eq.(JT3|) to be 
positive (or, in the case of complex uj, to have positive real part); the neg- 
ative cu-values are needed to construct real solutions. Then, the r + -modes 
describe vibrational motions turning counter-clockwise coupled to spin pre- 
cessions turning clockwise, i.e., opposite to the natural spin precession, and 
the r_-modes describe vibrational motions turning clockwise coupled to spin 
precessions turning counter-clockwise, i.e., in the same sense as the natural 
spin precession. 

When the lower sign is taken in Eqs.(|ll]), corresponding to a spin par- 
allel to the magnetic field, n = +z, we find two r + -modes with complex- 
conjugate mode frequencies. Thus, one of the mode frequencies possesses a 
positive imaginary part, indicating that the spin-up state is unstable for any 
value of K. We therefore concentrate on the stationary state with the spin 
antiparallel to the magnetic field, n = — z, corresponding to the upper sign 
in Eqs. (|TT|) . For the spin-down state we find one r + -mode for any value of 
K, and for K smaller than a critical value K c = a/4/27 two r_-modes, with 
real frequencies, as is shown in Fig.(^J). 

For K — > 0, the r + -mode, which goes like uj ~ uj v ^ (1 — K/2), and the 
slower of the r_-modes, which behaves as uj ~ uj vib (1 + K/2), become de- 
generate with frequency uj = uj vi b, corresponding to two linear independent 
purely vibrational modes. With increasing K, the coupling between trans- 
lational motion and spin motion lifts the degeneracy and gives rise to an 
increasing spin component, which leads to a decrease of the r + mode fre- 
quency and an increase of the slow T_ mode frequency. The fastest of the 
r_-modes, on the other hand, which for K —>■ is a pure spin precession 
mode with frequency uj = uj^/K = uj preci acquires with increasing K an in- 
creasing vibrational component, leading to a decrease of the mode frequency. 
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Figure 1: Real and imaginary parts of the mode frequencies as a function of 
K. 



At K = K c , it becomes degenerate with the slower r_-mode, with mode fre- 
quency lo = \/3u V ib. For K > K c , the two r_-modes have complex-conjugate 
frequencies. Fig. ([I]) shows the real and imaginary parts of the frequencies us 
of the three modes as a function of K. 

The dependence of the modes on the parameter K is shown more ex- 
plicitly by the form of the eigenvectors. The general form of these is given 
by 





( 


Bo ^vib 


\ 


( z ) - 




B\ id 








LO 




\ 


Ul v ib 


J 



where A± are dimensionless amplitude parameters. 



3.3 The excitation energy of the modes. 

The excitation energy of a given mode £ is defined as the difference between 
the energy of the mode and the energy of the stationary state, 



£ = —/in • B+-m 



d5x 
~~dt 



+ 



d5y 
~~dt 



(J,B . 



(21) 



Note that the energy contains bilinear terms in the coordinates and hence, 
one cannot neglect the ^-component of the spin. Instead, one must set 



n • z 



1 



(4 + 4) 



Thus, the correct expression of the energy for small amplitudes is 



- (e 2 x + e 2 ) B + B' ± (5xe x - 5ye % 



1 

-~m 
2 



d5x 



+ 



ddy 
~dT 



(22) 



The deviations 8x(t), Sy(t), e x (t), e y {t) from the stationary state belonging to 
the normal modes are given in real form by 



8x(t) 



-itat 



(>±m ( (■■<■■■ ty(t) = ±^p±, e~ - c.c. 



-iut 



(23) 
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,(t) = \t Tfi e- lut + ex., e y (t) = T Ie Ti0 e-^ + c.c. (24) 



With the help of Eq.(pD[), one obtains 



e = MBo ^ ^V *| 2 - (25) 
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ab 



From this result we conclude that for < K < a/4/27, the excitation energy 
of the vibrational modes, for which u 2 < Su>l ib , is positive while the excitation 
energy of the precessional mode, satisfying u 2 > 3u 2 ib , is always negative. 
At the point K = a/4/27, where the clockwise vibrational mode and the 
precessional mode coalesce, the excitation energy vanishes. We will further 
refer to these observations in the following section. 



3.4 The effect of viscous friction. 

When friction is introduced into the system, the equations of motion become 

d 2 x d dx , n . 

m ^ = / ^ (n ' B) - r ^ (26) 
d 2 V d dy 

m ^ = % {n - B) - n Tt 

and 

S— = an x B-r„nx — , (27) 
dt ^ p dt v 1 

where r t and r p are translational and precessional friction coefficients, respec- 
tively. The second term in the right-hand side of Eq. (|27j) is the spin-damping 
contributed by the change in the direction of the spin from n to n+rfn. Since, 
by definition, n is a unit vector, it must point perpendicular to dh. Thus, 
Q± = \dh/dt\ is the angular velocity associated with the change of n. Since 
the direction of Q± must be perpendicular to both dh and n we form the 
cross product Q± = fix (dh/dt) which incorporates both the correct value 
and the right direction. Multiplying Q± by r p yields the spin-damping term. 
To first order in r r and r t the secular equation in this case is given by 

= -K°u> + <4 - 2ul + 1 + 2iKojff - 2iKV/j (j^ * (28) 
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where we defined 



w n = 

to make the expression simple. Let u n fl be the eigenfrequencies u n of the 
frictionless problem, given by Eq. flT9|) . When adding small friction to the 
problem, the eigenfrequencies will change by a small amount 8u n . We find 
an approximate expression for Su n by expanding Eq. fl2"g| ) around u n to first 
order in 5u; n and making use of Eq. (p79|) . This gives 

^ \ W n,0 ~ 6 / 

Eq.(|2TS|) has an interesting consequence: The numerator in Eq.(|29[) is pos- 
itive for all three modes while the denominator is negative for the two vibra- 
tional modes and positive for the precessional mode. We therefore conclude 
that friction, either translational or precessional, stabilizes the vibrational 
modes and, simultaneously, destabilizes the precessional mode. The system 
all together becomes of course, unstable. 

The fact that spin damping leads to an exponential growth of the preces- 
sional mode is no surprise in view of its negative excitation energy. Also, the 
exponential decay of the vibrational modes due to translational friction is to 
be expected on account of their positive excitation energy. What is impor- 
tant is the fact that due to the coupling between translation and precession, 
translational friction causes an exponential growth of the precessional mode, 
with a growth time which, compared to the effect of spin damping, is smaller 
by a factor of r t K 2 S 2 / [imr v B^ in the limit of small K. 

4 Quantum-mechanical analysis. 

4.1 The Hamiltonian and its diagonalized form. 

In this section we consider the problem of a neutral particle with spin half ( 

5 = h/2) in a 2D inhomogeneous magnetic field from a quantum-mechanical 
point of view. Unlike the classical analysis, in which the derivation was valid 
for any value of the adiabaticity parameter K, we concentrate here on the 
behavior of the system when K is small. We choose to analyze the case 
of a spin half particle because this case already shows the essentials of the 
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quantum-mechanical problem. Note also that, quantum mechanically, the 
magnetic moment \i and the spin S of a particle are related by 



fi = jS, 

where 7 is the gyromagnetic ratio of the particle. Setting fi = 'jS and 
S = h/2 in Eq.(0) gives 



K 



'h(B 



I \2 



2 r )mB\ 



Now, it is convenient to express the spatial dependence of the magnetic 
field in polar coordinates r = y/ ' x 2 + y 2 and = arctan (y/x). We also 
denote by B the amplitude of B, by #-its direction with respect to the z axis 
and by <p the angle between the projection of B onto the x-y plane and the 
Thus, Eq.(0) is rewritten as 



B = B [sin 9 cos (p~x+ sin 9 sin ipy + cos 9z\ (30) 



where 



B\ x " 



B = B oX ll+[^) r 2 , (31) 



9 = arctan 



B'r 



B, 







{B v \ (y 
if = arctan I — - I = — arctan I — 



Thus, B and 9 depends only on r whereas (p depends (linearly) only on 
The time-independent Schrodinger equation for this system is 



h 2 

V 2 — iiB (sin 9 cos ipa x + sin 9 sin Lpa v + cos 9a z 

2m 



(r, 0) = EV' (r, 0) 
(32) 

where a x , a y and <j 2 are the Pauli matrices given by 

1 \ A / -z \ A f 1 

a «= l 1 ^ = U az= [o -1 
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E is the eigenenergy and is the two-components spinor 



ip\ (r, 0) 
ip\ (r, 4>) 



In matrix form Eq.(p2|) becomes 



I N i f , f ^'( r ^) 
^ii^ + -HmJ jit ,\ 



V4 ( r ' 

^1 0, 4>) 



where H K and if^-, given by 



h 2 2 

X 2m 



(33) 



(34) 



(35) 



M 



-liB 



cos 6* sin i 



-up 



sin i 



— cos 9 



are the kinetic part and the magnetic part of the Hamiltonian H, respectively. 

In order to diagonalize the magnetic part of the Hamiltonian, we make 
a local passive transformation of coordinates on the wavefunction such that 
the spinor is expressed in a new coordinate system whose z axis coincides 
with the direction of the magnetic field at the point (r, <p). We denote by 
R (r, <fi) the required transformation and set \l/ = Rty' . Thus, \l/ represent the 
same direction of the spin as before the transformation but using the new 
coordinate system. The Hamiltonian in this newly defined system is clearly 
given by RHR -1 . In the case of the magnetic field given in Eqs.([30l) and (|3l"l) 
the required transformation is accomplished by using the three Euler angles: 
First, we perform a rotation through an angle ip around the z axis. Second, 
we make a rotation through an angle 9 around the new position of the y axis. 
At the end of this process the new z axis coincide with the direction of the 
magnetic field. Now the value of the last Euler angle, which is a rotation 
around the new z axis, has no effect on this axis. For simplicity we choose 
this angle to be 0. Thus, the representation of the complete transformation 
for spin half particle is given by |1| 



R 



while its inverse is given by 





'.9, ' 








exp 


i~cr v 
. 2 . 


exp 


i—a z 
2 . 










.9 




exp 


—i—o z 
L 2 


exp 


—i- 
2 


fry 
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It is easily verified that the transformation indeed diagonalizes the magnetic 
part of the Hamiltonian as 



For the kinetic part we find, after some algebra, that 

1 ueV 1 

4 \ dr J 4r 2 

i „ „ d 



RHkR 1 — 

2m 



—i 



H — -a z cosQ 
r z a<p 

d6\ d_ 1<P6 ]_d9_ 

dr ) dr 2 dr 2 2r dr 

i - ■ a 9 
— -a r sin 6— - 



r* d(f) 

Thus, the Hamiltonian of the system in the rotated frame may be written 



as 



H = H diag + H, 



where 



Hdiag 
Hint — 



2m 
2m 



de 

dr 



l int 



1 i , n d 



Ar 2 



d<p 



/iBa z 



(36) 



(37) 



d£\ d_ + 1<P0_ l_d£ 
dr ) dr 2 dr 2 2r dr 



—<j x sm9— 

r z 0(p 



The first part of the Hamiltonian H^ag is diagonal. It contains the kinetic 
part ~ V 2 , a term whose form is =F fiB which is to be identified as the 
adiabatic effective potential and the terms ~ l/r 2 ,ir~ 2 a z d/d(f) which appear 
due to the rotation. The second part of the Hamiltonian H int contains only 
non-diagonal components. These will be shown to be of order O (K) and 
hence may be regarded as a small perturbation. We proceed to find the 
eigenstates of H diag . 
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4.2 Stationary states of H^iag- 

Since Hdiag is diagonal, the two spin states of the wavefunction are decoupled. 
We then seek a solution of the form 







; E = E 



(38) 



referred to as the spin-down state, and another solution 

1, = (Mr,*)) iE = Elt 



(39) 



which we call the spin-up state. 

The equation for the non-vanishing component of the spin-down state is 
given by 



2m 



v2 _ 1 fdtiy _ 1 _ i cosQ d 
4 \dr J 4r 2 r 2 86 



+ f iB\iJ l = E l iJ l , (40) 



whereas the equation for the non-vanishing component of the spin-up state 
is 



2m 



i fdey 1 z 8 

7 7 - — + -^COSV— 
4 V dr } 4r z r z am 



HB)^ = E^ V (41) 



We now show that in the limit of small K we can neglect the term ~ 
(dO/drf in both Eq.([p and Eq.(0): We compare the order of magnitude 
of the term fiB to that of the term h 2 (d6/dr) 2 /8m. Using Eq.(|3"TD it can 
be easily shown that the maximum value of d9/dr is B'j_/B whereas the 
minimum value of fiB is nB . Thus, 



h 2 fde 



lm V dr 



8{imB% 2 
(B' ± ) 2 h 2 = K 2 ' 



and hence we can neglect the term ~ (d6/dr) 2 when K is small. Furthermore, 
as we are interested in the solutions near the origin we replace the cos# 
term by its zeroth-order approximation around r = 0. We will justify this 
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approximation later. Under these approximations, Eqs. (|4*0| ) and (|¥ID simplify 

to 



2m 



v 2 1 * ^ 

4r 2 r 2 9(f) 



(42) 



and 



2m 



i d 



(43) 



The approximate solutions of these equations is outlined in the next two 
subsections. 



4.2.1 Stationary spin-down states. 



Eq. (|42j ) represents a particle in a symmetric attractive potential. If the extent 
of the wave function is small enough we can expand B in Eq.(pT|) to second 
order in r 



B~Bn 



1 (B\r 



+ O (r 4 ) 



(44) 



and apply the well-known solution of the harmonic oscillator 2(| in two di- 
mensions. Under this approximation, Eq. (|42"l) becomes 



h 2 



i d if 



2m \ r dr dr 2 r 2 



d 1 



V(B> ± ) 2 T 2 



+ 



2Bn 



We seek a solution whose form is 

^(r,0) = /(r)e^ 
and then the equation satisfied by / (r) is 



$1 = ( E l~ ^ B o)^i- 
(45) 

(46) 



2m 



r dr dr 2 r 2 



v — 



+ M ^ V/ = (^-MBo)/, (47) 
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which is an eigenvalue problem for /. The smallest eigenvalue for this prob- 
lem is obtained by setting 



v = — 



1 

2' 



for which the eigenfunction / is 



/ (r) = D exp 



imi (B' ± 
Ah 2 B 



D exp 



1 ( B\r 



4K V B 



Thus, under the harmonic oscillator approximation the down-part of the 
spin-down state is 



B\ 



B V2nK 



exp 



1 (B\r 



AK V B, 



Ad> 2 



(48) 



where the normalization constant D has been calculated by demanding that 

oo 2n 

rirjd* 



using the definite integral 



re ar dr 



2a 



Note that the extent of this wave function over which it changes appreciably 
is given by 



Ar, ~ \[K 



Bo 



B 



(49) 



whereas the extent over which fiB changes significantly (see Eq. (PT|) ) is 

Ar MB ~ (50) 
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Thus, the ratio between these two length scales is 



Ar, 



Ar 



K. 



(51) 



We therefore conclude that when K is small enough, the harmonic approx- 
imation is justified. Note also that Ar^ is also the typical length of cos 9. 
This shows that the substitution of cos 6* in Eq. (f40"D by 1 is also justified. 



The wave function ip^, given by Eq.(|8), then represents the lowest possi 



ble bound state for this system. This state corresponds to a trapped particle. 
The energy of this state is clearly 



El = A^o + 2 I —uj vi b 



fxB (l + 2K) ~ fiB 



(52) 



while its full spinor representation is 
/ 



\ 



\ B V2nK 



exp 



AK V B, 



B\r 



J<t>/2 



(53) 



4.2.2 Stationary spin-up states. 

Eq.(^) describes a particle in a repulsive potential. It corresponds to an 
unbounded state representing an untrapped particle. In this case there is a 
continuum of states, each with its own energy. As we are interested in non- 
radiative decay, we focus on finding a solution with an energy which is equal 
to the energy found for the trapped state, that is 

£ T = E t ~ nB . (54) 

When evaluating the lifetime in the next section, we compute the matrix 
element of H int between the states ip^ and ip^. Thus, most of the contribution 
to this integral comes from the region in r where ip^ is substantial. According 
to Eq. (|51|) , fiB changes very little in this range and, as a first approximation, 
we may take cos 9 ~ 1 and the potential in this region as uniform, 

fiB ~ fiB (55) 

in Eq. fl4"3]) . We now set a solution whose form is 

V> T (r,0)=s(r)e^. 
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Substituting this, together with Eqs.(Q) and flSSP into Eq.(^) gives 



2m 



Idg d 2 g g 



2/j,B g, 



whose non- singular solution is 



9{r) 



J. 



7+1/2 



B' ± r 



where J a (x) is the Bessel function of the first kind of order a. 

We note that all the four terms of Hi nt does not operate on the <fi co- 
ordinate. This is a consequence of the fact that L z — S z (where L z is the 
z-component of the orbital angular momentum and S z = ha z /2 is the z- 
component of the spin) is conserved. Hence, in order to have a non-vanishing 
matrix element between the up-state and the down-state, they must have the 
same 0-dependence. Thus, 7 = v — 1/2 and as a result 



"0| = CJ X 



with 



^1 



B' ± r 



B' ± r 

~B~oK 




J<t>/2 



(56) 



(57) 



where C is the normalization constant which is chosen to be real. 

The wave function given in Eq. fl5§) is oscillatory. It has a period of about 



near the origin. Comparing it to Arj_ given in Eq.(^9|), we find that 

Ar, 



(58) 



(59) 



which shows that, for K <C 1, the wavef unction ip^ executes many oscillations 
in the region where ipi is appreciable. 
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4.3 The lifetime. 



To evaluate the lifetime T esc of the particle in its trapped state, which is the 
average time it takes for the particle to escape, we calculate the transition 
rate from the bound state given by Eq. fl53|) , to the unbounded state Eq. fl5T|) , 
according to Fermi's golden rule^TJ. Thus, 



T 



2n \H I 2 



where 



2tt 



H u = JrdrJ d^\H int ^ j 







fde\ 


9 . 


ld 2 6 


1 d6 i 


[dr J 


dr 


2dr 2 


2r dr r 2 



sin 9 







(60) 



(61) 



is the matrix element of H int Eq.(|37|) between and and g(-&r) is the 
density of the final states at energy E^. 

The integrand in Eq.([H]) consists of a product of three elements: The 
function whose 'width' is about Ar^ (given in Eq. (|49|) ) around the origin, 
An operator consisting of four ^-dependent terms whose extent around the 
origin Ar M # is roughly \Jl/K larger than Arj_ and the function ip^ which 
is an oscillatory function with a characteristic period near the origin Arf 
which is smaller than Ar^. This suggests that we can approximate the 
integral in Eq. (|6l|) by substituting sin#, dd/dr and d 2 6/dr 2 by their value at 
r = 0, 



B' 

sin 9 ~ — — r 

B 

d9 B\ 
dr B 
d 2 9 

~ 

dr 2 

Substituting Eqs.(H), (HD and flU) int o Eq.(U) gives 



(62) 



H 



* a -^ h2 ^k c n exp 



1 

K 



(63) 
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where we have used the definite integral 



r 2 Ji (br) e 



dr= — exp 



4a 



(64) 



When Eq.fl63|) is substituted into Eq.© the term C 2 #(£ T ) appears. This 
term can be calculated by temporarily introducing suitable boundary con- 
ditions: Assume that the system is bounded by an infinite potential wall at 
r = R, the radius R being large compared to Ar j yet small when compared to 
Ar^B- In this case the uniform potential approximation holds for all r < R, 
and the wave function has the form 

^ (r, ( j>) = Cg(r)e^ 2 , 

where the radial part g(r) satisfies the Schrodinger equation 

h 2 



2m 



d 2 g 1 dg 
dr 2 r dr 



y 



(jiB + E)g. 



This equation has as non-singular solutions the Bessel functions of order 1, 

J\{kr) where k 2 



9{r) 



2 777, 



with eigenvalues k = k n determined by the boundary condition. From R 
rj it follows that kR 3> 1 such that J\(kR) may be approximated by the first 
term of its asymptotic expansion. Therefore, the boundary condition reads 



J^kR) 



7T 



kR 



cos I /, /)' — — ) 

4 / 



0. 



This yields the eigenvalues k n = (n + 1/4)tt/R. The density of states on the 
fc-axis is thus given by dN/dk = R/n, from which one obtains the density of 
states on the energy axis at E = fiB 



p{E = nB ) 



dN 
~dE 



mR 



1 

2^ 



rn 



R. 



nh 2 k 2tt Y h 2 fxB 
The constant C is determined by the normalization condition 



(65) 



\^\ 2 rdrd(j) = 2nC 2 / [J^kr^rdr 



R 
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which gives 



f R [J 1 {kr)} 2 rdr= l -R 2 [J 2 {kR)f. 
Jo * 



In the asymptotic region kR 3> 1, the function J 2 (kR) takes the values 
±^/2/(rrkR) at the zeros of J\{kR). This gives 



and therefore 



C ~2R 



C 2 p(E = f,B 



m 
2^2' 



Finally, using Eqs.(^) and (|63D inside Eq. (|60|) gives 



T 



647ru;„i fe 



exp 



2 

K 



T, 



i lb 



128tt : 



exp 



2 



(66) 



(67) 



where T vi b = 2ir/u) V ib is the period of classical oscillations inside the trap. 



5 Discussion. 

Summarizing all we have found we conclude that the problem we have studied 
has three important time scales: The shortest time scale is T prec , which is 
the time required for one precession of the spin around the axis of the local 
magnetic field. The intermediate time scale is T v n, = T prec /K, which is the 
time required to complete one cycle of the center of mass around the center 
of the trap. These two time scales appear both in the classical and the 
quantum-mechanical analysis. The longest time scale (provided K is small) 
T esc , which is not present at the classical problem, is the time it takes for the 
particle to escape from the trap. 

Whereas the classical analysis yields an upper bound of K = ^/4/Z7 for 
trapping to occur, no such sharp bound exists in the quantum-mechanical 
analysis. This is related to the fact that one cannot associate an effective 
potential well with a finite barrier with the system. 

As an example, we apply our results to the case of a neutron and an 
atom trapped with a field B = 100 Oe and B /B' ± = 10cm. These pa- 
rameters correspond to typical traps used in Bose-Einstein condensation 
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experiments f2"B], ^J. The results, being correct to within an order 

of magnitude, are outlined in the following table. 





B = 100 Oe 




B /B' ± 


= 10cm 




Neutron 


Atom 


m gr 


~ io- 25 


~ IO" 22 


(j, emu 


~ io- 23 


~ 10- 20 


K 


~ 10~ 5 


~ 10- 8 


Tprec SGC 


~ 10" 6 


~ 10- 9 


T v i}y SGC 


~ 10" 1 


~ 10" 1 




~ io( loB ) 


~ io( 108 ) 



We note that in both cases K is much smaller than 1. Also, the calculated 
lifetime of the particle in the trap is extremely large, suggesting that the 
particle (either neutron or atom) is tightly trapped in this field. 

The problem studied in this paper deals with a spin 1/2 particle. Though 
this fact has little influence on the solution of the classical problem, the 
extension to higher spin values complicates the analysis of the quantum- 
mechanical problem. In this case one has to deal with a (2S + l)-component 
spinor, and the interaction Hamiltonian does no longer connect the (— S)- 
state to the (+S')-state, but only to the (— S + 1) and (-S + 2) states which 
for S > 5/2 will still be trapped. 

References 

[1] A. L. Migdall, J. V. Prodan, W. D. Phillips, T. H. Bergeman and H. 
J. Metcalf, "First observation of magnetically trapped neutral atoms", 
Phys. Rev. Lett, 54 (24), 2596-2599 (1985). 

[2] T. Bergeman, G. Erez and H. J. Metcalf, "Magnetostatic trapping fields 
for neutral atoms", Phys. Rev. A., 35 (4), 1535-1546 (1987). 

[3] V. S. Bagnato, G. P. Lafyatis, A. G. Martin, E. L. Raab, R. N. Ahmad- 
Bitar and D. E. Pritchard, "Continous stopping and trapping of neutral 
atoms", Phys. Rev. Lett, 58 (21), 2194-2197 (1987). 



25 



[4] W. Petrich, M. H. Anderson, J. R. Ensher and E. A. Cornell, "Sta- 
ble, tightly confining magnetic trap for evaporative cooling of neutral 
atoms", Phys. Rev. Lett, 74 (17), 3352-3355 (1995). 

[5] M. O. Mewes, M. R. Andrews, N. J. Van-Druten, D. M. Kurn, D. S. 
Durfee, W. Ketterle, "Bose-Einstein condensation in a tightly confining 
DC magnetic trap", Phys. Rev. Lett, 77(3), 416-419 (1996). 

[6] The Levitron is available from 'Fascinations', 18964 Des Moines Way 
South, Seattle, WA 98148. 

[7] The U-CAS is available from Masudaya International Inc., 6-4, Kura- 
mae, 2-Chome, Taito-Ku, Tokyo, 111 Japan. 

[8] R. Harrigan, U.S. Patent Number: 4,382,245, Date of Patent: May 3, 
1983. 

[9] Hones et al, U.S. Patent Number: 5,404,062, Date of Patent: Apr. 4, 
1995. 

[10] R. Edge, "Levitation using only permanent magnets" , Phys. Teach. 33, 
252-253 (1995) and "Corrections to the levitation paper", ibid. 34, 329 
(1996). 

[11] M. V. Berry, Proc. R. Soc. Lond. A 452, 1207-1220 (1996). 

[12] S. Gov and S. Shtrikman, Proc. of the 19 th IEEE Conv. in Israel, 184- 
187 (1996). 

[13] M. D. Simon, L. O. Heflinger and S. L. Ridgway, Am. J. Phys. 65 (4), 
286-292 (1997). 

[14] S. Gov, S. Shtrikman and H. Thomas, Los-Alamos E-Print Archive, 
[http://xxx.lanl.gov/physics/9803020l (1998). 

[15] S. Earnshaw, Trans. Cambridge Philos. Soc. 7, 97-112 (1842). 

[16] The quantized motion of atoms in a quadrupole magnetic trap has been 
studied numerically by T. H. Bergeman, P. McNicholl, J. Kycia, H. 
Metcalf and N. L. Balazs, "Quantized motion of atoms in a quadrupole 
magnetostatic trap", J. Opt. Soc. Am. B, 6 (11), 2249 (1989). Here, we 
use an analytic method to find the lifetime of the particle in such a trap. 



26 



S. Gov, S. Shtrikman and H. Thomas, Los-Alamos E-Print Archive, 
[http: / / xxx. lanl.gov / physics / 9808007] , Submitted to Am. J. Phys. 

J. D. Weinstein and K. G. Libbrecht, Phys. Rev. A., 52 (5), 4004-4008 
(1995). 

L. D. Landau and E. M. Lifshitz, Quantum Mechanics (Pergamon 
Press), 3 rd ed., pp. 213-214. 

'Quantum Mechanics' by E. Merzbacher, John Wiley & Sons., 2 nd Ed., 
Ch. 5, Sec. 3, 57-61. 

'Quantum Mechanics' by E. Merzbacher, John Wiley & Sons., 2 nd Ed., 
Ch. 18, Sec. 8, 475-481. 

This integral may be found in 'Table of Integrals, Series, and Prod- 
ucts' by I. S. Gradshteyn and I. M. Ryzhik, Academic Press, 5 th Ed., 
6.521.1, pp. 697. Note that this integral make explicit use of the fact 
that Ji (kR) = 0. 

M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman and 
E. A. Cornell, "Observation of Bose-Einstein condensation in a dilute 
atomic vapor", Science 269, 198 (1995). 

K. B. Davis, M. O. Mewes, M. R. Andrews, N. J. van Druten, D. S. 
Durfee, D. M. Kurm and W. Ketterle, "Bose-Einstein condensation in 
a gas of Sodium atoms", Phys. Rev. Lett. 75, 3969 (1995). 

C. C. Bradley, C. A. Sackett, J. J. Tollett and R. G. Hulet, "Evidence 
of Bose-Einstein condensation in an atomic gas with attractive interac- 
tions", Phys. Rev. Lett. 75, 1687 (1995); ibid 79, 1170 (1997). 

E. A. Cornell and C. E. Wiemann, "The Bose-Einstein condensate" , Sci. 
Am. 278, 26-31 (1998). 



27 



